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Abstract 

An asymptotic expansion of the two-loop two-point "master" diagram with two masses 
m and M, on the mass shell Q 2 = M 2 , is presented. The treatment of the non- analytical 
terms arising in the expansion around the branching point is discussed. Some details of 
the calculation of a new class of two-loop integrals are given. 
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1 Introduction 



Phenomenological problems in particle physics often require calculations of multi-loop 
Feynman diagrams involving fields of various masses. Since exact results for diagrams 
with more than one mass scale practically do not exist beyond one loop, it is reasonable 
to apply some analytical approximations. In particular, one successfully applies explicit 
formulae for asymptotic expansions in various limits of momenta and masses ffl (see 
for an informal review). The underlying idea is to utilize a hierarchy of mass scales in the 
given problem to reduce it to a calculation of one-scale diagrams. There are two general 
classes of problems which can be solved using those expansions: the large momentum case 
(when some external momenta are much larger than other relevant mass scales) and the 
case of very heavy virtual particles. In both classes important phenomenological problems 
have recently been solved by asymptotic expansions . See also || where propagator and 
vertex two-loop diagrams were systematically calculated in various regions of momenta 
and masses. 

There is, however, an important class of diagrams with a heavy particle both inside and on 
its mass shell on the external legs, when off-shell formulae |IJ are generally not applicable^. 
Explicit formulae for asymptotic expansions in momenta and masses in some typical limits 
with large on-shell momenta were presented in || . In our paper we discuss the two-point 
on-shell integrals with two mass scales shown in Fig. 1 (m < M). 




(a) (b) 

Figure 1: Examples of two-point on-shell (Q 2 = M 2 ) diagrams with two mass scales; solid 
lines denote a heavy particle with mass M, and dashed lines - a light one with mass m. 

The one-loop diagram in Fig. la is a good example which illustrates both the difficulties 
of the previous approaches and the idea of the new method. Of course the exact result is 
well known: 



(k 2 — m 2 )[(k + Q) 2 — M 2 } 




in D/2 ( - + 2 — In M 2 — x 2 \nx — 2aV4 - x 2 arctan J^—- I W 



3 Although these formulae are guaranteed at least off the mass shell, they are also valid in some on-shell 
situations, for example, in the pure large mass limit. 
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with x = m/M. In this paper we use dimensional regularization with D = 4 — 2e; it is 
understood that the masses in propagators have an infinitesimal negative imaginary part 
which is not displayed. We drop Euler's 7 in the formulae. 

Suppose now we want to find the expansion of G\ around x = without calculating the 
above two-scale integral. Unfortunately it is impossible to apply just the Taylor expansion 
of the integrand around m = 0; the expanded propagator fc2 _ m2 = -^J2^=o [ w) l ea ds 
to infrared divergences. This difficulty is connected with the impossibility of the naive 
change of order of expansion and integration. General formulae of asymptotic expansions 
give, for this diagram, two contributions: besides this 'naive' Taylor expansion in m one 
has a contribution which is the Taylor expansion of another factor of the integrand, k i+ 2 kQ > 
around fc 2 = ||. (In the off-shell situation, this extra contribution would be [J] just 
the Taylor expansion of ^ + q- ) 1 2 _ M 2 around k — 0, which happens to be impossible now 
because Q is on-shell, Q 2 = M 2 .) 

We therefore expand around k 2 = and obtain = ^ £~ (-^)". The 

resulting integrals are homogeneous in Q and depend on one mass scale m only. It is easy 
to convince oneself that in the result the divergences generated in the Taylor expansions 
are canceled and the sum of both contributions leads to 



G\ = in 



D/2 



I 7T 

- + 2 - In M 2 - nx + x 2 (1 - In x) + -x 3 + 0(x 4 ) 

2n+l 



(2) 



We notice that the terms tt Ij^j appear even though m is found in the integrand 
only in an even power. These non-analytical terms arise because we expand the diagram 
around the point m = which is also a branching point (threshold) of this Feynman 
diagram. 

For the two-loop diagram in Fig. lb no exact formula is known. In this paper we give 
its asymptotic expansion in the limit m/M — > 0. The next section is rather formal; it 
describes the general formalism of the expansion. In section |3| this machinery is applied 
to the diagram lb; it turns out that there are four contributions (rather than two, as in 
the one-loop case). The calculation of a new type of integrals which arises in this context 
is discussed in some detail. 



2 Explicit formula of asymptotic expansion 

We consider the Feynman integral corresponding to a graph T when the masses Mj and 
external momenta Qi are large with the respect to small masses rrii and external momenta 
qi. Let us suppose that the momenta are non-exceptional. Let the large external momenta 
be on the mass shell, Q 2 = M 2 . Then the asymptotic expansion in the limit Qi, Mi — ■> 00 
takes the following explicit form ||: 

F T (Q h M h q u mi]e) ^ ^M 1 F T (Q h M h qi, e). (3) 

7 
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Here the sum is over subgraphs 7 of T such that 

(a) in 7 there is a path between any pair of external vertices associated with the large 
external momenta Qf, 

(b) 7 contains all the lines with the large masses; 

(c) every connectivity component jj of the graph 7 obtained from 7 by collapsing all the 
external vertices with the large external momenta to a point is 1PI with respect to the 
lines with the small masses. 

In general 7 can be disconnected. One can distinguish the connectivity component 70 
which contains the external vertices with the large momenta. 

The operator that is involved in the sum is a product Oi -M 7i of operators of Taylor 
expansion in certain momenta and masses. These operators are by definition applied to 
integrands of Feynman integrals over loop momenta. For connectivity components 7, other 
than 70 the corresponding operator performs Taylor expansion of the Feynman integral 
_F 7i in its small masses and external momenta. (Note that its small external momenta 
are generally not only the small external momenta of the original Feynman integral but 
also some loop momenta of T.) Consider now Ai 10 . The component 70 can be naturally 
represented as a union of its 1PI components and cut lines (after a cut line is removed 
the subgraph becomes disconnected; here they are of course lines with the large masses). 
By definition M. lQ is again factorized and the Taylor expansion of the 1PI components of 
70 is performed as in the case of other connectivity components 7,. 

It suffices now to describe the action of the operator M. on the cut lines. Let I be such a 
line, with a large mass Mj, and let its momentum be Pi+h where Pi is a linear combination 
of the large external momenta and k\ is a linear combination of the loop momenta and 
small external momenta. If P\ = Qi then the operator Ai for this component of 7 is 



By T x we denote the operator of the Taylor expansion in x around x — 0. In all other 
cases, e.g. when P t = 0, or it is a sum of two or more external momenta, the operator M. 
reduces to the ordinary Taylor expansion in small (with respect to this line considered as 
a subgraph) external momenta, i.e. 

1 1 

\ h + Pl) 2_ M 2 = - {Kk + Pl) 2_ M 2 

In all cases apart from the cut lines with Pf = Mf the action of the corresponding operator 
«M is graphically described by contraction of the corresponding subgraph to a point and 
insertion of the resulting polynomial into the reduced vertex of the reduced graph. 



(4) 

,=] 



(5) 



3 



3 The two-loop master diagram 



Let r be the two-loop self-energy graph (see Fig. lb) containing three heavy lines with 
the mass M, two light lines with the mass m and let the external momentum be on the 
mass shell, Q 2 = M 2 : 

G = rr f^fi m 

2 J J (k 2 - m 2 )(P - m 2 )(k 2 + 2Qk) [(k + I) 2 + 2Q(k + I)} (I 2 + 2QI) ' 1 ; 

According to the general formula (|3|) three types of subgraphs contribute to the asymptotic 
expansion in the limit m/M — > 0: 

(i) the graph T itself; 

(ii) the subgraph 71 consisting of the left triangle and the third heavy line as well as 
symmetric subgraph 72; 

(hi) the subgraph 73 consisting of three heavy lines. 
For (i) we expand the integrand around m = 0: 

d D kd D l 1 

(k 2 + 2Qk){{k + I) 2 + 2Q{k + l))(l 2 + 2QI) m " {k 2 - m 2 )(l 2 - m 2 ) 

00 n 

= E^E^n (7) 

n=0 j=0 

where 

rr d D kd D l 

} ~ J J (p)«+i(/2)6+i(p + 2Qk) [(k + I) 2 + 2Q{k + I)} (I 2 + 2QI) ' ( ' 

The leading term (which gives G2 at m = 0) is equal to [[7] 



7(0,0) 



7T 



D r g 



M 2 



|C(3)-vr 2 ln2 



+ 0(e). (9) 



For any a and b the integral 7(a, b) can be reduced to 7(0, 0) and one-loop integrals using 
recurrence relations || M. This algorithm has been implemented in FORM 

Each of the two contributions to (ii) equals 



d D l _ 1 r d D k _ 1 



T l 2 n , o^, / 13. , *^lJ rn*T3 -3 T l ,, , A9 , , A ( 10 ) 



I 2 - m 2 1 I 2 + 2QI J k 2 + 2Qk m k 2 -m 2 (k + I) 2 + 2Q(k + I) 



.73 



£ (_l)Ji m 2«i+iO £ 

jlJ2,j3=0 n=0 



'J3 

v n , 



r D (2Ql+m 2 ys-n , (2fcQ" 
7 {I 2 - m 2 )(2Ql)n+ l J (F)^+i(A; 2 + 2QA:y3+2' 1 J 
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The factors I 2 in the numerator have been substituted by m 2 because whenever I 2 — m 2 
appears it cancels the corresponding factor in the denominator and we get a massless 
tadpole integral which vanishes. The calculation of both one-loop integrals in (|Tl"|) poses 
no difficulties. 

Finally, the contribution (iii) is 

d D kd D l _ 1 



(k 2 - m 2 )(P - w?) ( K k 2 + 2Qk)(n(k + I) 2 + 2Q(k + l))(nl 2 + 2QI) 



K=l 



= y [f d n kdDl (-m 2 )n+n Hk + l)T 

h ,^ 3= JJ (k 2 - m 2 )(2Qk)^+ 1 (2Q(k + l))^(l 2 - m 2 )(2QZ) j2+1 ' 1 ! 

In this formula the following new type of integrals arises: 

J(oi, 02,03, 0.4? 05) = 

d D kd D l (13) 

(k 2 -m 2 + iO) a i {I 2 -m 2 + i0) a * (2Qk + i0) a ^ (2Qk + 2QI + i0) a * ' v ; 

In general we can also have a power of 2QI as an extra term in the denominator but this 
can be removed by partial fraction decomposition. Any such integral can be calculated 
analytically. We are interested in the case 01 = a-i — 1. In calculating fllcf) we first reduce 
the problem to the case 05 = by expressing the product {kl) ar ° in terms of traceless 
products (kip 1 ' = k^ a ' l 'lr a iy Then we notice that the result of the integration over I of a 
traceless product l( a ,i) times the part of the integrand which depends only on / must be 
proportional to the traceless product Q( a ,i)- Therefore the factor (kl)^' can be replaced by 
(Qkp l >(Qiy % ' I (QQY l \ Then the factors involved are expressed through ordinary products 
Qk and Ql. 

Now, an arbitrary integral J(ai, 02, 03, 04, 0) is reduced to integrals with 03 = with the 
help of the following relation obtained by integration by parts 

(2D - 2 - 2a 4 - a 3 )3- - (D - 1 - a 4 )4~ + 4a 3 3+ = 0. (14) 

In the present problem we only need the integrals J with the unit values of the first two 
arguments Oi = a 2 = 1. In general it is also easy to reduce a 12 to 1 using a similar 
recurrence relation. 

The relation (|14]) can be used to express an arbitrary integral J(l, 1, a 3 , a 4 , 0) through 
integrals with a 3 = 1 and a 3 = 0. The latter will be considered below. To calculate 
J(l,l,l,a4,0) we use the symmetry with respect to I <-> k and replace 2QI in the nu- 
merator by Qk + Ql, which decreases a 4 . In the resulting integral we use partial fraction 
decomposition with respect to 2QI. As a result we get integrals with decreased a 4 and 
products of one-loop integrals. If 03 becomes increased in this process we apply the 
relation (111) again. 



5 



In the result of manipulations described above, the only two-loop integrals we are left 
with are of the form J(l, 1, 0, a, 0) = J (a) 



J(a) J J d D kd D l^ 2 _ m2 + iQ ^ p _ m2 + iQ ^ 2 Qk + 2QI + ' ^ 

In J(a) we first perform the ko and 1$ integrations using Cauchy theorem; the angular 
integrations are trivial since there are no products kl in the integrand. The two remaining 
radial integrations lead to 



C 2 2 x ~ 2e coslvre + f /3 a \ (a a 3 \, , 

J (a) = —. — V . 2 J B (-_ e -_l + e ]B(--2 + 2e, — + - - e J 16 

v ; (-4Mm) a cos(tt£) V2 ' 2 ) \2 ' 2 2 ) y ' 

with C = 2m ^3 g?r ^ — . This formula is valid to all orders in e. 

4 The final result 

After adding the contributions (i), (ii), and (hi) we obtain the following finite result 



4 r o 18 ~| 

G 2 = — |7r 2 ln2 - |C(3) + £ [a« + + ln(x 2 )c n ] x n + 0(x 19 )\ . (17) 



The coefficients a n and c n are given in Table 1. The computing time for these first 18 
powers of x is of the order of a few hours on a DEC ALPHA workstation. For the 
coefficients of 7r 2 we find &2n+i = 0, &2„ = — If this is true for all n this part of the 
series can be summed up analytically: tc 2 J2%Li &2n^ 2n = tt 2 ln(l — x 2 /2). 

Similarly to the one-loop result @ we see in the expansion ( |P7D non-analytical terms with 

/ \2n+l / \ 

odd powers of x; at two-loop we also obtain logarithmic terms 7r ( || J In f ^j. 

The expansion (|T^) converges at least up to m/M = 1. We have compared the values 
obtained with this formula with the numerical program For x < 0.5 the accuracy 



is better than 10 . Even at x — 1 the error is only 0.7%; at this point we can use for 



comparison the recently found analytical value of the curly bracket in eq. (|T7|) at x — 1 



-C(3) + |vrCl2(7r/3) ||12|| . Much more important than these numerical achievements is the 



possibility of extending the expansion (T7\) to any order in x. The methods described in 
this paper can also be extended to other cases of propagators and to vertex diagrams. For 
example, the results of ref. for two-loop QCD corrections to b — > c transitions at zero 
recoil can be extended to other kinematical regions. Work on this is in progress. 



J We thank J. Franzkowski for sending us an independent evaluation of G2 at x = 1. 



6 



n 


a n 


c n 


1 





71 


2 


-1/2 





3 


2tt/3 


tt/24 


4 


-5/144 


-1/12 


5 


13tt/60 


3tt/640 


6 


107/5400 


-1/15 


7 


5297r/6720 


5tt/7168 


8 


9073/470400 


-73/1680 


9 


14887tt/483840 


35tt/294912 


10 


101923/7938000 


-17/630 


11 


715801tt/56770560 


63tt/2883584 


12 


3620783/461039040 


-1391/83160 


13 


31515089tt/5904138240 


231tt/54525952 


14 


47658179/10100170080 


-317/30030 


15 


3278369671vr/1416993177600 


143tt/167772160 


16 


24219791/8480609280 


-19741/2882880 


17 


13114766971tt/12847404810240 


6435tt/36507222016 


18 


8503364437/4825201661280 


-21071/4594590 



Table 1: Coefficients of the expansion of the diagram G2 
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